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Abstract 

An almost model-independent parametrization for the ratio of the total cross section to the 
elastic slope, as function of the center of mass energy, is introduced. The analytical result is based 
on the approximate relation of this quantity with the ratio R of the elastic to total cross section 



and empirical fits to the R data from proton-proton scattering above 10 GeV, under the conditions 
7^ , ' of asymptotic unitarity and the black-disk limit. This parametrization may be useful in studies of 

extensive air showers and the determination of the proton-proton total cross section from proton-air 
production cross section in cosmic-ray experiments. 
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I. INTRODUCTION 



In addition to their intrinsic astrophysical importance, cosmic-ray experiments constitute 
a valuable tool for the investigation of particle and nuclear physics at energies far beyond 
those obtained in accelerator machines. However, at the highest energies a direct approach 
to particles properties and their interactions is difficult due to the decreasing flux with the 
increase of the energy. Presently, an indirect method, based on extensive air shower (EAS) 
studies, is the usual way to treat the subject [l]. In these events, the distribution of the 
first interaction point allows, in principle, the determination of the proton-air production 
cross section ^ and in a second step, the estimation of the most fundamental quantity in 
hadronic interactions: the proton-proton total cross section 0, 4|. 

However, in practice, the interpretation of the EAS development depends on extrapola- 
tions from phenomenological models that have been tested only in the accelerator energy 
region, resulting therefore in systematic theoretical uncertainties. That represents a cru- 
cial point because different models with distinct physical pictures and at the same time 
consistent with the experimental data up to cm. energies ~ 2 TeV, present, in general, 
contrasting extrapolations at the cosmic-ray region (above ~ 50 TeV). Therefore the theo- 
retical uncertainties (bands) involved are relatively large and very difficult to be estimated. 
As a particular consequence, the estimations of the proton-proton total cross section from 
different cosmic-ray experiments and analyses are characterized by large error bars and even 
discrepant results, as discussed in js] and references therein. 

EAS studies are essentially based on Glauber's multiple diffraction theory 6|, |7| and its 
extensions and/or corrections, including Gribov-Regge screening effects and other ingredi- 
ents. In this context the evaluation of the hadron-nucleus elastic and quasi-elastic cross 
sections, which contribute to the nucleon-air production cross section, depends on the ratio 
of two physical quantities, the total cross section atot and the elastic slope B. That ratio 
just represents one of the main sources of uncertainties in model extrapolations. 

In this work an almost model-independent parametrization for the ratio (Jtot/B is pro- 
posed, which may avoid uncertainties from models tested only at lower energies. The 
parametrization is based on the approximate, but experimentally justified, connection of 
(Jtot/B with the ratio R of elastic to total cross section. By means of unitarity arguments 
and empirical fits to R data from pp scattering above 10 GeV and up to 7 TeV, an analytical 
parametrization for R{s) is introduced and then extended to the ratio atot/B. 



The manuscript is organized as follows. In Sect. II we introduce the physical quanti- 
ties of interest with explicit reference to the importance of the ratio (Jtot/B in cosmic-ray 
studies. In Sect. Ill we recall some formal (rigorous) results from axiomatic QFT and how 
some inequalities can be connected with experimental results. In Sect. IV we introduce 
the analytical parametrization and present the fit results. The conclusions and some final 
remarks are the contents of Sect. V. 

II. PHYSICAL QUANTITIES AND THE GLAUBER FORMALISM 



Let us first recall the main physical quantities related to high-energy elastic hadron scat- 
tering, defining the notation and normalizations jsj. Neglecting spin effects and denoting 
F(s, t) the invariant elastic amplitude in terms of the Mandelstam variables, s and t, the dif- 
ferential and total cross sections at high energies (s >> 1 GeV^) are expressed, respectively, 
by 



da , , IGtt , ^ , , ,9 
-(s,t) = -^\F(s,t)\\ 



(1) 



and 
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ImF(s,t = 0) (Optical Theorem). 



(2) 



The parameter p, the ratio between the real and imaginary parts of the forward amplitude, 
is given by 

ReF(s, t = 0) 



Pis) 



ImF(s, t = 0)' 

and the slope of the elastic differential cross section in the forward direction is defined as 



(3) 



B{s,t = 0) 



d / da\ 
di \~dt) 



t=o 



From (1) to (3), the optical point is expressed by 



da 
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The integrated elastic cross section reads 



° da 
, dt^ 



s, t)dt, 



(4) 



(5) 



(6) 



where to defines the physical (kinematic) region and, from unitarity, the inelastic cross 
section is obtained: 

O-m(s) = (Ttot{s) - ael{s). (7) 

For our purposes let us recall two formulas in the Glauber formalism that play a central 
role in EAS studies The first one is the expression for the sum of the elastic and 
quasi-elastic cross section for hadron-nucleus (hA) scattering, 




where fj and bj are the coordinate and impact parameter of the individual nucleons, T{fj) 
the single nucleon density, 6 the impact parameter of the cosmic-ray hadron and aj{s, b — bj) 
the nucleon-nucleon impact parameter amplitude (profile function). In addition to possible 
configurations for the nucleus, the profile function constitutes the main ingredient for the 
connection between hadron-hadron and hadron-nucleus scattering. Typically this profile is 
parametrized by [l| 



a,(.,6,) = i^^:^^e-^l/[^^W] 



An 



B s) 



(9) 



where p, atot and B are the quantities defined above, demanding inputs from models to 
complete the connection. As clearly illustrated by Ulrich et al. the uncertainty bands 
for these three quantities resulting from high-energy extrapolations based on representative 
phenomenological models, are larger than the range covered by all available MC interaction 
model results, as QGSJETOlc, EP0S1.61, S1BYLL2.1 and QGSJETII.3, leading the authors 
to the conclusion that presently, "the extrapolation of hadronic cross sections to cosmic-ray 
energies might be underestimated" 

In what concerns the above three fundamental quantities, we recall that forward am- 
plitude analyses connect atot{.s) and p(s) through dispersion relations. Detailed tests on 
different parametrizations have been developed by the COMPETE Collaboration, with the 



selection of the highest rank result 
Physics by the Particle Data Group 



10| . which also appears in the Review of Particle 



Recent results by the TOTEM Collaboration on 



atot at 7 TeV [12| and the expected estimation of p at this energy, will certainly shed light on 



novel analytical parametrizations and therefore more reliable extrapolations with less model 
dependency. 

However, that is not the case for the elastic slope B{s) since the available data from 
pp and pp elastic scattering can be extrapolated in a very large band of possibilities and 
moreover, any result is strongly model dependent. As a consequence, from Eq. (9), despite 
the uncertainty in the effective radius of the nucleon-nucleon amplitude, any extrapolation 
is strongly dependent on the ratio 

^{s), (10) 

namely the unknown correlation between atot and B in terms of energy. Although some 
analytical connections have already been investigated from fits to the experimental data 



13| . the statistical and systematic errors in both quantities and the model dependencies 
involved put limits on these results. 

Based on the above comments, we understand that even under some reasonable approx- 
imate conditions, an almost model-independent parametrization for the above ratio may 
reduce the uncertainty band in the extrapolations from accelerator to cosmic-ray-energy 
regions. That is the point we are interested in here. 



III. FORMAL AND EXPERIMENTAL RESULTS 

In this section we first recall some rigorous (formal) results related to crtot{s) and B{s) 
and their connections with the experimental data presently available. Based on these con- 
siderations in the next section we introduce our proposed parametrization and present the 
results. 



III.l. Rigorous Results 

General principles and high-energy theorems have always been a fruitful source of model- 
independent results for physical quantities in the asymptotic regime [l^. In this context, 
two well-know inequalities have been derived for the total cross section and the elastic slope. 
The first one is the Froissart-Lukaszuk-Martin upper bound, stating that asymptotically 
(s — )■ oo) 
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for some Sq ISNlTj. The second one, playing here an important role, is the lower bound of 
MacDowell and Martin, obtained from unitarity together with properties of the Legendre 
polynomials and involving forward quantities [ISj, 



^lnImF(s,t) 



> 



1 <tis 



(12) 



j^^g ISvr aei{s) 

From the definition of the forward slope, Eq. (4), together with Eq. (1) and under the 
assumption 



it follows that 



ReF{s,t = 0) « lmF{s,t = 0), 



dt 



In ImF(s,t) 



-B 



t=o 



and from (12) an upper bound is obtained for our ratio of interest. 



— — < IStt — . (13) 

B{S) (Ttotis) 

Although rigorous, relations involving inequalities have a limited practical applicability, 
except for bounds imposed on the construction of phenomenological models. To go further in 
the search for almost empirical or model-independent results, experimental data and formal 
inequalities must be checked, as follows. 

III. 2. Experimental Results 

The highest energies reached in accelerator machines for particle and antiparticles reac- 
tions concern pp and pp scattering, covering the region up to ~ 2 TeV (pp) and, presently, 
up to 7 TeV (pp). These data indicate that at the highest energies 

pis)<O.U, 



which means that the above assumption, p << 1, constitutes a considerable approximation. 
On the other hand, at the optical point, Eq. (5), an assumption like 

1 + p2 ^ 1 (14) 

certainly represents a reasonable approximation. We shall return to this point in what 
follows. It may be interesting to note that these information allow us to derive bound (13) 
from (12) under different assumptions, as shown in Appendix A. 

Concerning the differential cross section, experimental data indicate a sharp forward 
peak, followed by a dip-bump or dip-shoulder structure above ~ 0.5 GeV^ (Tevatron, LHC). 
Typically, these structures are located more than 5 decades below the optical point, Eq. (5). 
These experimental facts are important in the determination of the integrated elastic cross 
section, since in this case the differential cross section can effectively be represented by an 
exponential fall off, simulated by a model-independent parametrization jl2|. 



da da 
dt da"^ 



(15) 



with B the (constant) forward slope. In that case, with the reasonable approximation (14) 
at the optical point (5) and assuming to — t- — oo in Eq. (6), the integrated elastic cross 
section reads 



O-elis) 



B{s) 167r 

and therefore, 



ajM ^ (16) 

B{s) atot{sy ^ ' 

which is very close to the approximate bound (13). However, the main ingredient in this 
result is the possibility to investigate the behavior of atot{s)/B{s) from formal and experi- 
mental information on the ratio aei{s)/atot{s), as discussed in what follows. 

IV. ANALYTICAL PARAMETRIZATION AND FIT RESULTS 

In Figure 1 we display the experimental information presently available on the ratio 



o^ei/o'tot from pp scattering above 10 GeV [ll|, including the recent TOTEM result at 7 TeV 
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12| (highest energy reached in accelerators). From a strictly empirical point of view, the 
data in the linear-log scale may suggest a parabolic parametrization in terms of In s, with 
positive curvature. However, unitarity demands an obvious bound, 

— < 1. 

In addition, naive models, as for a Gaussian profile or the gray-disk, predict 13] (Jei/atot = 
C/2, where C is a constant (absorption coefficient) and in the black-disk limit C = 1. These 
results indicate a constant asymptotic limit for the ratio 

lim = A (constant) 

and therefore a change of sign in the curvature, at some finite value of the energy, is expected. 
Moreover, since from Fig. 1 the data at low energies indicate (Jei/cTtot ~ constant ^ 0.18, a 
general behavior related to a logistic or sigmoid function can be conjectured, at least above 
10 GeV (the high-energy region). Several functions with this property can be considered. 
However, for tests on goodness of fit some quantitative information on the value of the 
asymptotic limit A is necessary. 

Looking for a wide range of possibilities in the phenomeno logical context, we shall consider 



two contrasting pictures that have been discussed in the literature. On t 
amplitude analysis by Block and Halzen favour the asymptotic black-disk 



le one hand, the 



20|, which is also 



predicted, for example, in the models by Chou and Yarrg M arrd by Bourrely Softer and 

Wu |22]. On the other hand, the U-matrix unitarization scheme by Troshin and Tyurin 
predicts atotis) ~ CTeiis) ~ In^ s and crj„e/(s) ~ Ins ^], which is beyond the black-disk limit 
and in agreement with the above obvious unitarity bound. Therefore these two contrasting 
pictures suggest 

A = - (black-disk limit) and - < A < 1 (beyond the black-disk limit). 

Although, in principle, it might be possible to explore all the real interval for A beyond the 
black-disk, we consider here only its maximum value. As we shall show, that is adequate 
and sufficient for our purpose to infer wider uncertainty bands from the above mentioned 
phenomenological context. 

Based on these considerations and inspired in different physical phenomena, we have 
tested several functional forms to fit the aei/<Jtot data. The best statistical result has been 
obtained with the following novel model-independent parametrization: 

8 



— (s) = Atanh(7i + 72 Ins + 73 111^3), (17) 

(7 tot 

where 7^, z = 1, 2, 3 are free fit parameters and A represents the asymptotic limit, for which 
we consider only the two extreme cases A = 1/2 and A = l. 

The data reductions have been performed with the objects of the class TMinuit of ROOT 
Framework [24]. We have adopted a Confidence Level of ~ 68 % (one standard deviation), 
which means that the projection of the distribution in [N + l)-dimensional space {N = 
number of free fit parameters) contains 68 % of probability 25|. The fit results for A = 1/2 
and A = 1 are displayed in Fig. 2 and Table I together with the statistical information. In 
both cases, the error propagation from the fit parameters has been evaluated and are also 
displayed in the figure; the bands however are indistinguishable. 

From the approximate result (16), the ratio crtot/B can be predicted as function of the 
energ y an d in an almost model-independent way. The results together with the experimental 



data 



11 



12 



26| are displayed in Fig. 3 and show that, in fact, Eq. (16) is very close to 



the approximate bound (13). Up to our knowledge, the only rigorous result indicating a 




Vs (GeV) 



FIG. 1. Experimental data on the ratio between the elastic and total cross sections from pp 
scattering above 10 GeV |ll|,Q. 
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constant asymptotic value for this ratio appears in the recent formal analysis by Azimov, 
on boundary values for the physical cross section and slope j2^. The numerical predictions 
with uncertainties for the ratio (Jtot/B at the LHC energy region are displayed in Table II, 
together with the experimental value at 7 TeV. 

TABLE I. Fit results with parametrization (17) for the ratio CTei/cTtot from pp scattering above 10 
GeV. In both cases the degrees of freedom (DOF) are 87. 





A = 1/2 


A = 1 


7i(xlO-i) 


4.66± 0.18 2.204± 0.078 


72(xlO-2) 


-2.59±0.49 


-1.11±0.20 


73(xlO-3) 


1.77±0.33 


0.76 ±0.13 


xVdof 


1.167 


1.168 




10 10^ 10^ 10" ^10^ 10^ 10^ 10^ 10^ 

Vs (GeV) 



FIG. 2. Ratio between the elastic and total cross section and fit results through parametrization 
(17), including uncertainty from error propagation, for A = 1 (upper curve) and ^ = 1/2 (lower 
curve) . 
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TABLE II. Predictions from Eqs. (16-17) for the ratio atot/B at the LHC energy region and the 
TOTEM result at 7 TeV nJ- 





A = 1/2 


A = 1 


TOTEM 


7.0 TeV 


12.827±0.047 


12.821±0.024 


12.56±0.59 


14 TeV 


13.811±0.068 


13.903±0.033 





V. CONCLUSIONS AND FINAL REMARKS 



Based on unitarity arguments and fits to the experimental data on the ratio R = crf.i/atot 
from pp scattering above 10 GeV, a novel empirical parametrization for R{s) has been 
introduced, Eq. (17). The approximate connection between this quantity and the ratio 
atot/B, Eq. (16), allows us to infer the corresponding energy dependence for this ratio in an 
almost model-independent way. All the results are in agreement with rigorous inequalities 




?1G. 3. Experimental data on the ratio between the total cross section and the elastic slope 
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12 



26l | and predictions from Eqs. (16-17), including uncertainty from error propagation, for 



A = 1 (upper curve) and ^ = 1/2 (lower curve). 
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derived from the axiomatic formulation. 

Although depending on the unknown asymptotic limit represented by the constant A, 
the results here presented lead to, at least, four main conclusions: 



1. If the black-disk represents a reliable physical limit 20|, its saturation is very far from 



presently available energies: ^/s > 10^ GeV, from Figure 2; 

2. If the Froissart-Lukaszuk-Martin bound is saturated then in this region B(s) ~ In^ s; 

3. Either for A = 1/2 or A = 1, the uncertainty bands evaluated by error propagation 
from the fit parameters in Eq. (17) are negligible: upper, central and lower curves in 
Figures 2, 3 and 4 overlap; 

4. Even with A = 1/2 and A = 1 as lower and upper bounds, extrapolation of the ratio 
(Jtot/B to the Auger energy region, ^/i ~ 50 - 60 TeV, indicates a reasonably small 
error band, overestimated from Figure 4 to be in the interval 15.5 — 16.3. 

At last, we understand that the applicability of our results in the context of the Glauber 
formalism can be further developed and improved along the following lines: 




FIG. 4. Detail of Fig. 3 up to the Auger energy region. 
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The recent TOTEM result 
lined in Subsection III. 2 



"or aei at 7 TeV has been obtained through the steps out- 
12| and therefore has been evaluated from the results for atot 
and B. The forthcoming measurement of atot, by means of a luminosity-independent 
method, and the corresponding CTg/ determination may improve our fit result, since 
this region is just associated with the change of curvature in parametrization (17), as 
shown in Fig. 2. Moreover, further results at 14 TeV will certainly contribute with 
additional improvements in the fit results. 

Here we have limited the discussion to the extreme values indicated from unitarity and 
the black-disk limit, together with the almost model-independent result for the ratio 
<^tot/B. However, as commented in Sect. HI, analytical parametrizations for the total 
cross section and the p parameter, as those obtained by the COMPETE Collaboration 



(or possible deviation from this result, if confirmed by the experimental data 28|). 
may be combined with our results in the Glauber connection, Eqs. (8-9), reducing the 
uncertainties bands in the extrapolations to cosmic-ray energies. 
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Note added 

After this paper was submitted for publication, we have noticed the results from a gray- 
disk- model analysis, in which a different transition on o'ei/o'tot{s), from low to high energies, 
is proposed ^29| (see also references therein). 
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Appendix A 

Beyond the forward direction we can define 



, da , , 



and 



p(s,t) 



RoF(,s.O 
ImF(s,i)' 



so that from Eq. (1), 



t) = 2^ lnImF(s, + | ln[l + p\s, t)\. 



Under the reasonable assumption that at least in the neighborhood of t=0 

ImF(s, t) > ReF{s,t) 



and by expanding the second term in the r.h.s of (A2) we obtain at t—0 



t=o 



+ o{pHs)). 



t=o 



Since from the experimental data p{s) < 0.14 and under the assumption 

d 



Eq. (Al) at i = reads 



lim-p(.,t) = 0, 



B(s) w 2— lnImF(s,t) 
^ ^ dt ^ ^ 



t=o 



leading through Eq. (12) to the upper bound (13). 
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